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Abstract 

In this article, the estimation of rehabihty of a system is discussed p{y < a;)when 
strength, X, and stress, Y , are two independent exponential distribution with differ- 
ent scale parameters when the available data are type II Censored sample. Different 
methods for estimating the reliability are applied . The point estimators obtained 
are maximum likelihood estimator, uniformly minimum variance unbiased estimator, 
and Bayesian estimators based on conjugate and non informative prior distributions. A 
comparison of the estimates obtained is performed. Interval estimators of the reliability 
are also discussed. 
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1 Introduction 



In life testing, it is often the case that items drawn from a population are put on 
test and their times of failure are recorded. For a number of reasons, such as budget or 
limited time, it is often necessary to terminate the test before all the failure times have 
been observed; In this case the data become available in some ordered manner. The test 
is usually terminated after a fixed time or a fixed number of failures are observed giving a 
censored sample see Epstein and Sobel (1953). In stress-strength model, the stress Y and 
the strength X are treated as random variables and the reliability of a component during 
a given period (0,T)is taken to be the probability that the strength exceeds the stress 
during the entire interval, the reliability of a component is P{Y < X). Several authors 
have considered different studies for stress and strength with exponential distribution with 
complete sample. Tong(1974) discussed the estimation ofP{Y < X) in the exponential 
case. Tong(1977) had a look at the estimation of P(Y < X) for exponential families. A 
good review of the literature can be found in Johnson (1988). Beg(1980a,b and c)estimated 
the exponential family ,two parameter exponential distribution and truncation parameter 
distributions . Basu(1981) considered maximum likelihood estimators (MLE) for P{Y < X) 
in case of gamma and exponential distributions. Sathe and Shah (1981) studied estimation 
of P{Y > X) for the exponential distribution. Chao (1982) provided simple approximations 
for bias and mean square error of the maximum likelihood estimators of reliability when 
stress and strength are independent exponentially distributed random variables. Awad and 
Charraf (1986) studied three different estimators for reliability has a bivariate exponential 
distribution. Dinh, et al (1991) obtained the MVUE of R when X and Y have the bivariate 
normal distribution. Moreover, they considered the case when X and Y have the bivariate 
exponential distribution. Bai and Hong (1992) estimated P{Y < X) in the exponential case 
with common location parameter. Kunchur and Mounoli (1993) obtained UMVUE of stress- 
strength model for multi component survival model based on exponential distribution for 
parallel system. Siu-Keung Tse and Geoffrey Tso (1996) studied the shrinkage estimation 
of reliability for exponential distributed lifetime. A note on the UMVUE on P(Y < X) in 
the exponential case discussed in Cramer and Kamps (1997). Selvavel, et al (2000) studied 
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reliability R when (X, Y) jointly follows a truncated bivariatc exponential distribution with 
a common parameters. Khayar (2001) discussed the reliability of time dependent stress- 
strength models for exponential and Rayleigh distributions. Shrinkage estimation of P{Y < 
X) in the exponential case discussed by Ayman and Walid (2003).Tachen (2005)deals with 
the empirical Bayes testing the reliability of an exponential distribution . 

In the present article, the reliability, R , is studied when X and Y two independent 
exponential distribution with different scale parameters. Which can be represented as 
f{x;a) = , x > 0,Q! > 0. 

Different estimators of R are derived, namely, maximum likelihood estimator (MLE), uni- 
form minimum variance unbiased estimator, (UMVUE), and Bayesian estimators with mean 
square error loss functions corresponding to conjugate and non informative priors. A com- 
prehensive comparison of the various point estimators (MLE,UMVUE, and Bayes) is per- 
formed on the basis of the mean squared error. Interval estimators of R are also discussed. 
A numerical comparison of the intervals obtained. 

2 Reliability 

Let X be the strength of a component and Y be the stress acting on it. Let X and Y be 

exponential independent random variables with parameters a and /?, respectively. That is 
, the probability density functions (pdfs) of X and Y are, respectively, 

1 -X 

f{x; a) = -e « , x > 0, a > (2.1), 

and 

fiy;P) = ^e^, y>o,p>o (2.2), 

where a and /? are unknown parameters . 
The reliability of the component will be 

R = P{Y < X) 
= ^ / Ooo(l - e^)e^dx 

= fracaa + (i. (2.3) 
if a and (3 are known then R is simply calculated using Eq.(2.3). 
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3. Point Estimation of R 



3.1 Mciximum Likelihood Estimator of R 

If a and /? are unknown the MLE of, i?i , of i? is given by 

Ri = fracaa + /3, (3-1) 

where a and (3 are the MLEs of a and (3 , respectively. For obtaining a and we argue as 
follows: 

Suppose that ri components where ri < n with strengths Xi;i = 1, ...,ri , each of which 
having exponential distribution with parameter a as in Eq. (2.1) are subjected, respectively. 



the likelihood function will be 

(n-ri) 



77,1 1 - l^i-l ^'i -^ri (ri-ri) 

L = —^i-re^r' h ,egno(3.2) 

(n — rij! a 



Taking the logarithm of Eq. (3.2) and find the derivative with respect to a 

-na = Xi + {n- ri)xri , 

^^E£i^. + (n-n)x,,^ (3.3) 

To stress Y}; j = 1, r2, having exponential distribution with parameter (3 asm. Eq. (2.2), 
where r2 < m . Assuming that Xi and Ij; z = 1, ...,ri and j = 1, ...,r2 , are independent , 



the likelihood function will be 

m\ TT 
(m-r2)! -l^i 



E^2 
i^ = r^^(ire^ (3.4) 

(m — r2)! p 



Taking the logarithm of Eq. (3.4) and find the derivative with respect to/? 

(UnL _ -r2 , X]j=i I {m-r2)yr2 _ n 
~ - ~ + ^ 0^ - ^ 

-r2P = E7=i + (rn - r2)yr2, The MLE of /3 will be 



^^J-Uyj + im-r,)y.,^ (3.5) 
r2 



Now we shall study some properties of Ri . 
First, if ri = r2 = r : 

1] Unbiasedness 

^(i^i) = ^[1-^1^(1-^)^] 

' r+1 r+1 

limr_>oo E{Ri) = R - i^linir^oo 
then lim^^oo = R 

then, i?i asymptotically unbiased estimator of R . 
2] Consistency 

Var{R,) = ^^[^n^^]' 

lim^^oo 1^or(i?i) = R^^ lim^^oo f 
then 

\im.y^^Var{Ri) = 

then, R\ is a consistent estimator for R. 

Second, if ri r2 ■ 

1] Unbiasedness 

For fixed r2, 

lim,,_oo£^(i2i) = lim,,^oo ^^^[1 - Klir7-2) i^ - 
lim,,_oo£^(i2i) = R[l - ^(1 - i?)'] 
then limr^^r2^c» 

then, Ri asymptotically unbiased estimator of R . 
2] Consistency 

Vnr(R.\ — (^1+^2-1) r ('■2-1)" I2r 1 12 

('•2-l)/3 ^+(ri-l)a 

For fixed r2. 
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lim^j^oo Var{Ri) = lim^i_oo rlt?-2) ^ ^ij^^m^? limri->oo[— ]^ 

and, 



limri,r2-»oo Var{Ri) = limrj-^oo ^ 

then 

yar(i?i) = 
then, Ri is a consistent estimator for R. 

3.2 Uniform Minimum Variance Unbiased Estimator of R 

Let Xi, Xr^SLudYi, ...,1^2 be two independent random samples , of size ri andr2 , respec- 
tively, drawn from exponential distributions with parameters a and P , respectively. 
Define 

Zi = lne^',Vj = Ine^^ ,i = l,...,riandj = l,...,r2, 
Z = J2i = lri2;i,andF = J2j = 1^2 

Clearly from Eq. (3.2) and (3.4) we see that Z,V is a complete sufficient statistic for a , p. 
Now, we have 

E{W) = l.P{vi < Zl) + O.P{vi > Zl) 
E{W) = P{lneyi < Ine^i) 
E{W) = P{y < x) = R 

LetT^ be the indicator variable /[q 2^)(Ty).it could be seen that W be unbiased estimator 
for R , by using Rao-Black Well and Lehmann-Scheffe we have R2 is UMVUE for R .(see 
Mood et al (1974)). 
R2 = E{W/Z, V) 

^2 = wf{zi,vi/Z, V)dvidzi 

where f{zi,vi/Z,V) is the conditional pdf of -21,^1 given Z,V . Notice that ziand vi are 
independent exponential random variables with parameters a and /?, respectively, and that 
Z and V are independent gamma random variables with parameters (n, a) and (m, (3) , 
respectively. 

Wc sec that Z — Zl and V — vi are impendent gamma random with parameters (n — l,a) 
and {m — 1,P) , respectively. Moreover Z — zi and zi ,as well as V — vi and vi are also 
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independent. Wc see that 

Rn-f r ,,. r(ri)(Z-^i)('-l-^)r(r2)(V-»;i)('-2-2) 

' Tin - l)r(r2 - I m^-^{Z-z,)(^^-^)(y-vi)(^^-^Uv,dzu v, > zv, 

(3.6) 

The computation of the UMVUE R2 is very comphcated as it can seen from equation 
(3.6) .so, we will use the MATHCAD program to evaluate the value of R2 . 

3.3. Bayes Estimator of R 

We obtain Bayes Estimator of R with respect to the mean square error loss function with 
respect to conjugate and non informative prior distributions. 

3.3.1. Conjugate gamma prior distribution 

Let Xi, ...jXr^andYi, ...jYj.^ be the first ri and r2 failure observations from Xi, and 
Yi, ...jYm respectively, where both of them have exponential distribution with parameters 
aand/3 respectively. Assume that the prior distribution of a is given by 

vroi = /(a) = ^^{^)"'-^e^,ui,vi,a> 

the likelihood function with type II censored sample is , respectively, 

f{x,,...,XrAa) = (l)ne-^(E.=i'-ia.,+..,(n-n))^ (3.7) 

[n — ri)\ a 



and 



f{yi,...,yM = ^^(i)'-^e-ME,.,-^.+^^.(— )), (3.8) 



Assuming that aand/3 are independent having prior gamma distributions, the posterior 
distributions of aand/3 will be gamma distribution also, 

I , (^1 + E£i x^ + xr, (n - ri))"i+-i (i)'-^+"^e- . 
-.2 = f{a\x„ = r(n + m + ' 

(3.9) 
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and 

7i"3 = /(/3|yi,-,yr2) 



r(r2 + n2 + 1) 
the joint posterior function, put, 



(3.10) 



ri r2 



i=l j=l 
^ = r(ri+«i+l)r(r2+M2+l) 

7r(a,/?|x,y) = A;(("i+'-i)r("2+'-2)(i)n+«i^ ^r2+«2g^g^^ (311) 

a p 

Hence Bayes estimator of R will be 

l./-(ui+ri) {u2+r2) nl f?{ri+ui)(-i _ f>\{r2+U2+l) 

Rs = EiR\x, y) = ^, , : , ^, , / tT^^^tA^tStT^TT^^^^ (3-12) 



r(ri + ui + r2 + U2 + 1) Jo ((1 - ii;)C + Et)«i+''i+«2+'-2 

Prom equation (3.12) there is no explicit form of i?3SO, The computation of the Bayes 
estimator i?3,is very complicated ,we will use the MATHCAD program to evaluate the 
value of i?3 . 

3.3.2. Non Informative Prior Distributions 

Let Xi, ...,Xri be a random sample from exponential distribution with parameter a . The 



prior distribution of a is proportional to y/I{a) , where I{a) is Fisher's information of the 
sample about a , and is given by 



/(a) = ^, (3.13) 



from that the prior distribution 



TTl oc -, (3.14) 

a 



Similarly, if Yi, Y,../is a random sample from exponential distribution with parameter p, 
the prior distribution of /3 will be given by: 

772 oc ^, (3.15) 
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if we have a and /3 are independent then the posterior joint distribution of a and /3 ,will be 



n{a,(3\xi,...,Xri,yi,...,yr2) L{xi, |/3)7ri(a)7r2(/3), (3.16) 

then 

7ria,P\x,,...,Xr„yi,...,yr,) = , ^, (i)-ie-^g,=inx,+x.,(n-n)). 



{n — T\)\ a {m — r2)\ 



put 



<^ = ^^rxXi^Xr^{n - ri)),e = (^r2yj + yriim - r2)) 

7r(a,/:/a;i,...,Xri,yi,...,yr2) = 7 ttt >«>P > 0> (3-17) 

(n — ri)!(m — r2)' 

under the mean square error, Bayes estimator R4 of R will be 



n\m\ 



1 i2(^2+2)(i_i2)(n+i) 



i?4 = = — / i——j.dR, (3.18) 

^ ' '^^ (n-ri)!(m-r2)!(ri +r2 + 3)! 70 {{I - R)5 + Rey^+'-^+^ ' ^ ^ 

Prom equation (3.18) there is no explicit form of i?4SO, The computation of the Bayes 
estimator i?4,is very complicated ,we will use the MATHCAD program to evaluate the 
value of i?4 . 



4. Interval Estimation of R 

4.1 Approximate confidence interval 

Let Xi, Xj.^ and li, Y^.^ be a random samples with size ri, r2 from exponential distri- 
bution with parameter a, /? respectively ,we can show that the maximum likelihood function 
of R ,R is asymptotically normal distribution with mean R and variance-covariance matrices 
I~^{a,l3) where I{a,l3) is the Fisher information matrix and given by: 




9 



The variance-covariance matrix is obtained by inverting the information matrix with 
elements that are negatives of the expected values of the second order derivatives of log- 
arithms of the likelihood functions, and the asymptotic variance-covariance matrix is ob- 
tained by replacing values by their maximum likelihood estimators. Hence, the asymptotic 
variance-covariance matrix will be 

Prom (Surles and Padjett (2001)) that the maximum likelihood function of R ,Ri is asymp- 
totically normal distribution with mean R and variance 

2 rir2ia + Pf 
o', — 



Hence, (1 — q;)100 an approximate confidence interval for R would be (Li, Ui) ,and the value 
of Li, Ui is given as: 

Li = i?i-zi_|a^^, (4.1) 

and 

C/i = i?i + Zi_|a^^, (4.2) 
where quantile of the standard normal distribution and Ri is given by Eq.(3.1). 

4.2 Exact confidence interval 

Let Xi, ...,Xrj and Yi,...,Yr2 be random samples with size ri,r2 from exponential distri- 
bution with parameter a, [3 respectively. 

Since Zi = /ne^^and Vj = Ine^^ are independent with gamma distribution with parame- 

^^^^ y. ^, ^^^^ — i=i^' and 2(3lne^i=^^^ are independent with chi 

square distributions with degree of freedom 2ri and 2r2 hence rRi = {1 + ^)~^ 

we know that Fi = fracVaZ/3 has P-distribution with (2ri, 2r2)degrees of freedom ,rRi = 

(1 -I- ^Fi)~^ ,this equation can be written as follow Fi = ^'^^J^''^ (i-R) 

using Fi as a pivotal quantity, we obtain a (1 — a) 100% confidence interval for R as 

L2 = Fi_|(2r2,2ri)(Fi_|(2r2,2ri) + ^)-^ 
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and 

f/2 = F|(2r2,2ri)(F|(2r2,2ri) + 

where 

Fojis (1 — a)th of an F distribution random variables with (2r2, 2ri) degrees of freedom. 

5. Numerical illustrations 

in this section we will compare the different point estimators of R , namely Ri,R2,R3 and R4 and 
in the cases of R = 0.25, R = 0.4, R = 0.5 and R = 0.538 ,and different values of their pa- 
rameters a ,(5 . 2999 samples are generated of various size of n = 5, 10, 15, 20 , n = 25 and 
n = 50 from exponential distribution with parameters a ,j3 respectively. 

5.1. Numerical illustrations in the case ri = r2 in tables(l-4) 

1- MLE has the smallest mean square errors expect at some points of the UMVUE is the 

smallest, and for a = 2, /3 = 6 .both Bayes estimator and Non-Informative estimator are 
the smallest. 

2- For some values of n, m, ri,r2,Q; and (5 Bayesian estimator has advantage over UMVUE 
comparing the mean square errors, such that: 

* In table (1) in case a = 2, /? = 3 at n = m = 5, ri = r2 = 3. 

* In table (2) in case a = 2, /? = 6 for all values of n, m except at n = m = 5, ri = r2 = 4. 

* In table (4) in case a = (3 = 7 for the values of n = m = 15, ri = r2 = 12 . 

3- For all values of n, m, ri and r2 Bayesian estimator equal to Non-Informative estimator. 

5.2. Numerical illustrations in the case ri 7^ r2 in tables(5-7) 

1- MLE has the smallest mean square errors expect at some points of the UMVUE is the 
smallest .and for a = 2, /3 = 6 .both Bayes estimator and Non-Informative estimator are 
the smallest. 

2- For some values of n, m, ri, r2, a and f3 UMVUE has advantage over MLE comparing the 
mean square errors, such that: 
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*In table (5) in case a = 2, /3 = 3 for the values of n = 5, m = 4, ri = 3, r2 = 2 and in the 
case n = 10, m = 5, ri = 6, r2 = 5, also in the case n = 25, m = 10, ri = 20, r2 = 9 . 
*In table (6) in case a = 2, P = 6 the UMVUE is the smallest in the all cases expect for the 
values of n = 10, m = 25, ri = 10, r2 = 20 and in the case n = 15, m = 25, ri = 10, r2 = 20 
,also in the case and n = 50, m = 25, ri = 4, r2 = 24. 

*In table (7) in case a = (3 = 7 the UMVUE is the smallest in the all cases expect for the 
values of n = 5, m = 4, ri = 3, r2 = 2 and n = 5, m = 5, ri = 4, r2 = 3. 

3- For some values of n, m,ri,r2,a and /3 Bayesian estimator is equal to UMVUE comparing 
the mean square errors. 

4- all mean square errors (MSEl, MSE2, MSE3 gamma and MSE4) increases as a or /3 
increases 

when we compare the case of ri = r2 and ri ^ r2 we have that, in case ri ^ r2 the mean 
square errors of all points estimators are smaller than the mean square errors in r\ = r2 
case, expect some points. 



12 



Comparisons between different estimators for R = P{Y < X) in censored 
exponential CEise 

Table{l) 



m 


n 


n 


r2 


Ri 


MSEl 


R2 


MSE2 


i?3 


MSE3 




MSE4 


5 


5 


3 


3 


0.276 


0.0150 


0.536 


1.39 


0.00177 


0.159 


0.0000432 


0.160 


5 


5 


4 


4 


0.631 


0.0530 


1.90 


0.702 


0.00 


0.160 


0.00 


0.160 


10 


10 


6 


6 


0.62 


0.05 


0.66 


0.25 


0.00 


0.16 


0.00 


0.16 


10 


10 


7 


7 


0.70 


0.09 


2.14 


0.00 


0.00 


0.16 


0.00 


0.16 


10 


10 


8 


8 


0.44 


0.00 


0.58 


0.02 


0.00 


0.16 


0.00 


0.16 


10 


10 


9 


9 


0.76 


0.13 


0.63 


0.68 


0.00 


0.16 


0.00 


0.16 


15 


15 


12 


12 


0.70 


0.09 


0.73 


0.59 


0.00 


0.16 


0.00 


0.16 


15 


15 


13 


13 


0.59 


0.04 


0.58 


0.52 


0.00 


0.16 


0.00 


O.IG 


15 


15 


14 


14 


0.53 


0.02 


0.59 


0.57 


0.00 


0.16 


0.00 


0.16 


20 


20 


15 


15 


0.66 


0.07 


0.56 


0.53 


0.00 


0.16 


0.00 


0.16 


20 


20 


16 


16 


0.54 


0.02 


0.64 


0.75 


0.00 


0.16 


0.00 


0.16 


20 


20 


17 


17 


0.73 


0.11 


0.70 


0.00 


0.00 


0.16 


0.00 


0.16 


25 


25 


23 


23 


0.60 


0.04 


0.00 


0.61 


0.00 


0.16 


0.00 


0.16 


25 


25 


24 


24 


0.48 


0.01 


0.67 


0.73 


0.00 


0.16 


0.00 


0.16 


50 


50 


4 


4 


0.60 


0.04 


0.73 


0.73 


0.00 


0.16 


0.00 


0.16 


50 


50 


6 


6 


0.37 


0.00 


0.68 


0.76 


0.00 


0.16 


0.00 


0.16 


50 


50 


9 


9 


0.67 


0.08 


0.61 


0.45 


0.00 


0.16 


0.00 


0.16 



a = 2, (3 = SandR = 0.4 
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Table{2) 



ni 


11 


'■i 


'2 


-ffi 


MSEl 


R2 


MSE2 




MSE3 


Ri 


MSE4 


5 


5 


3 


3 


0.725 


0.226 


0.526 


0.076 


0.0003457 


0.062 


0.051 


0.04 


5 


5 


4 


4 


0.789 


0.291 


2.244 


3.978 


8.594E-11 


0.062 


0.0006678 


0.062 


10 


10 


6 


6 


0.773 


0.273 


0.668 


0.175 





0.063 





0.063 


10 


10 


7 


7 


0.657 


0.165 


0.789 


0.29 





0.063 





0.063 


10 


10 


8 


8 


0.671 


0.177 


0.656 


0.164 





0.063 





0.063 


10 


;io 


9 


9 


0.708 


0.209 


0.619 


0.1 36 





0.063 





0.063 


15 


15 


12 


12 


0.657 


0.165 


0.645 


0.156 





0.063 





0.063 


15 


15 


13 


13 


0.741 


0.241 


0.783 


0.284 





0.063 





0.063 


15 


15 


14 


14 


0.812 


0.316 


0.673 


0.179 





0.063 





0.063 


20 


20 


15 


15 


0.662 


0.17 


0.779 


0.28 





0.063 





0.063 


20 


20 


16 


16 


0.71 


0.212 


0.679 


0.184 





0.063 





0.063 


20 


20 


17 


17 


0.727 


0.228 


0.719 


0.22 





0.063 





0.063 


25 


25 


23 


23 


0.799 


0.301 


0.734 


0.234 





0.063 





0.063 


25 


25 


24 


24 


0.643 


0.154 


0.78 


0.281 





0.063 





0.063 


50 


50 


4 


4 


0.518 


0.072 


0.771 


0.272 





0.063 





0.063 


50 


50 


6 


6 


0.688 


0.192 


0.571 


0.103 





0.063 





0.063 


50 


50 


9 


9 


0.794 


0.296 


0.734 


0.234 





0.063 





0.063 



a = 2, (3 = 6andR = 0.25 
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Table{3) 



ni 


11 


'■i 


'2 


-ffi 


MSEl 


R2 


MSE2 


By, 


MSE3 




MSE4 


5 


5 


3 


3 


0.486 


0.002721 


0.896 


0.128 


0.0029 


0.287 


0.000016 


0.29 


5 


5 


4 


4 


0.546 


0.00006 


0.91 


0.138 


6.109E-08 


0.29 


1.652 


1.239 


10 


10 


6 


6 


0.516 


0.0005 


0.532 


0.000045 





0.29 





0.29 


10 


10 


7 


7 


0.607 


0.0046 


0.639 


0.01 





0.29 





0.29 


10 


10 


8 


8 


0.563 


0.00059 


0.569 


0.00095 





0.29 





0.29 


10 


;io 


9 


9 


0.55G 


0.00031 


1.703 


1.357 





0.29 





0.29 


15 


15 


12 


12 


0.434 


0.011 


0.578 


0.0016 





0.29 





0.29 


15 


15 


13 


13 


0.353 


0.035 


0.545 


0.00004 





0.29 





0.29 


15 


15 


14 


14 


0.413 


0.016 


0.616 


0.006 





0.29 





0.29 


20 


20 


15 


15 


0.429 


0.012 


0.532 


0.0000417 





0.29 





0.29 


20 


20 


16 


16 


0.543 


0.000023 


0.00045 


0.289 





0.29 





0.29 


20 


20 


17 


17 


0.482 


0.0032 


0.639 


0.01 





0.29 





0.29 


25 


25 


23 


23 


0.526 


0.00016 


0.0000053 


0.29 





0.29 





0.29 


25 


25 


24 


24 


0.46 


0.006215 


4.86E-12 


0.29 





0.29 





0.29 


50 


50 


4 


4 


0.484 


0.002912 


0.663 


0.016 





0.29 





0.29 


50 


50 


6 


6 


0.443 


0.009077 


0.72 


0.033 





0.29 





0.29 


50 


50 


9 


9 


0.453 


0.007321 


0.412 


0.016 





0.29 





0.29 



a = 7, 13 = QandR = 0.538 
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Table{A) 



ni 


11 


'■i 


'2 


-ffi 


MSEl 


i?2 


MSE2 


By, 


MSE3 


Ri 


MSE4 


5 


5 


3 


3 


0.535 


0.001238 


0.549 


0.002414 


0.007157 


0.243 


0.000089 


0.25 


5 


5 


4 


4 


0.286 


0.046 


0.592 


0.008417 


0.29 


0.044 


0.0000024 


0.25 


10 


10 


6 


6 


0.681 


0.033 


0.562 


0.003845 


0.00085 


0.249 





0.25 


10 


10 


7 


7 


0.519 


0.00035 


0.7 


0.04 





0.25 





0.25 


10 


10 


8 


8 


0.426 


0.00546 


0.713 


0.046 





0.25 





0.25 


10 


;io 


9 


9 


0.-112 


0.007753 


0.G28 


0.017 





0.25 





0.25 


15 


15 


12 


12 


0.525 


0.0006029 


1.619 


1.252 





0.25 





0.25 


15 


15 


13 


13 


0.628 


0.016 


0.055 


0.198 





0.25 





0.25 


15 


15 


14 


14 


0.682 


0.033 


0.6 


0.01 





0.25 





0.25 


20 


20 


15 


15 


0.43 


0.00485 


0.687 


0.035 





0.25 





0.25 


20 


20 


16 


16 


0.487 


0.0001711 


0.004066 


0.246 





0.25 





0.25 


20 


20 


17 


17 


0.426 


0.00552 


0.000001057 


0.25 





0.25 





0.25 


25 


25 


23 


23 


0.366 


0.018 


0.591 


0.008315 





0.25 





0.25 


25 


25 


24 


24 


0.46 


0.006215 


8.178E-14 


0.25 





0.25 





0.25 


50 


50 


4 


4 


0.516 


0.0002437 


0.647 


0.022 





0.25 





0.25 


50 


50 


6 


6 


0.487 


0.0001691 


0.659 


0.025 





0.25 





0.25 


50 


50 


9 


9 


0.443 


0.003232 


0.206 


0.086 





0.25 





0.25 



a = 7,P = landR = 0.5 
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Table{5) 



ni 


11 


'■i 


'2 


-ffi 


MSEl 


i?2 


MSE2 


Ry, 


MSE3 




MSE4 


5 


4 


3 


2 


0.677 


0.077 


0.362 


0.001465 


8.027E-09 


0.16 


0.00002938 


0.16 


5 


5 


4 


3 


0.521 


0.015 


0.759 


0.129 


2.095E-11 


0.16 


6.352E-08 


0.16 


10 


5 


6 


5 


0.78 


0.144 


0.457 


0.003292 





0.16 





0.16 


10 


10 


7 


6 


0.56 


0.026 


0.716 


0.1 





0.16 





0.16 


10 


15 


8 


7 


0.638 


0.057 


0.62 


0.049 





0.16 





0.16 


10 


20 


9 


10 


O.Gil 


0.011 


0.0001796 


O.IG 





O.IG 





O.IG 


10 


25 


10 


20 


0.524 


0.015 


0.000001812 


0.16 





0.16 





0.16 


15 


5 


6 


5 


0.796 


0.157 


0.567 


0.028 





0.16 





0.16 


15 


10 


7 


6 


0.533 


0.018 


0.8 


0.16 





0.16 





0.16 


15 


15 


8 


7 


0.546 


0.021 


0.621 


0.049 





0.16 





0.16 


15 


20 


9 


10 


0.518 


0.014 


0.601 


0.04 





0.16 





0.16 


15 


25 


10 


20 


0.589 


0.036 


0.001212 


0.159 





0.16 





0.16 


20 


5 


15 


4 


0.325 


0.005696 


0.65 


0.063 





0.16 





0.16 


20 


10 


16 


9 


0.614 


0.046 


0.59 


0.036 





0.16 





0.16 


20 


15 


17 


14 


0.551 


0.023 


0.531 


0.017 





0.16 





0.16 


25 


10 


20 


9 


0.724 


0.105 


0.689 


0.083 





0.16 





0.16 


25 


15 


22 


14 


0.657 


0.066 


0.644 


0.06 





0.16 





0.16 


25 


20 


24 


19 


0.569 


0.029 


0.574 


0.03 





0.16 





0.16 


50 


25 


-1 


24 


0.80G 


0.165 


0.226 


0.03 





O.IG 





O.IG 


50 


30 


6 


29 


0.727 


0.107 


0.018 


0.146 





0.16 





0.16 


50 


40 


9 


39 


0.759 


0.129 


0.0001917 


0.16 





0.16 





0.16 



a = 2,P = SandR = 0.4 
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Table{6) 



ni 


11 


'■i 


'2 


-ffi 


MSEl 


R2 


MSE2 


Ry, 


MSE3 




MSE4 


5 


4 


3 


2 


0.911 


0.436 


0.907 


0.432 


0.0003433 


0.062 


0.001571 


0.062 


5 


5 


4 


3 


0.715 


0.216 


0.811 


0.314 


1.118E-13 


0.062 


5.941E-07 


0.062 


10 


5 


6 


5 


0.805 


0.308 


0.875 


0.39 





0.062 





0.063 


10 


10 


7 


6 


0.837 


0.344 


0.825 


0.33 





0.062 





0.063 


10 


15 


8 


7 


0.727 


0.228 


0.739 


0.24 





0.062 





0.063 


10 


20 


9 


10 


0.71 


0.212 


0.73 


0.231 





0.0G2 





0.0G3 


10 


25 


10 


20 


0.762 


0.262 


0.0001145 


0.062 





0.062 





0.063 


15 


5 


6 


5 


0.865 


0.378 


0.789 


0.291 





0.062 





0.063 


15 


10 


7 


6 


0.702 


0.204 


0.8 


0.302 





0.062 





0.063 


15 


15 


8 


7 


0.801 


0.304 


0.59 


0.116 





0.062 





0.063 


15 


20 


9 


10 


0.792 


0.294 


0.776 


0.277 





0.062 





0.063 


15 


25 


10 


20 


0.729 


0.23 


0.003433 


0.061 





0.062 





0.063 


20 


5 


15 


4 


0.583 


0.111 


0.718 


0.219 





0.062 





0.063 


20 


10 


16 


9 


0.587 


0.114 


0.72 


0.221 





0.062 





0.063 


20 


15 


17 


14 


0.781 


0.282 


0.862 


0.374 





0.062 





0.063 


25 


10 


20 


9 


0.738 


0.238 


0.758 


0.258 





0.062 





0.063 


25 


15 


22 


14 


0.71 


0.211 


0.663 


0.17 





0.062 





0.063 


25 


20 


24 


19 


0.78 


0.281 


0.81 


0.314 





0.062 





0.063 


50 


25 


-1 


24 


0.841 


0.35 


0.395 


0.021 





0.0G2 





0.063 


50 


30 


6 


29 


0.818 


0.323 


0.889 


0.408 





0.062 





0.063 


50 


40 


9 


39 


0.859 


0.371 


0.00000394 


0.062 





0.062 





0.063 



a = 2,P = 6andR = 0.25 
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Table{7) 



ni 


11 


'■i 


'2 


-ffi 


MSEl 


R2 


MSE2 




MSE3 




MSE4 


5 


4 


3 


2 


0.61 


0.012 


0.436 


0.004035 


0.013 


0.237 


0.022 


0.229 


5 


5 


4 


3 


0.259 


0.058 


0.509 


0.00007698 


0.000102 


0.25 


0.009586 


0.241 


10 


5 


6 


5 


0.427 


0.005336 


0.602 


0.01 


0.0006316 


0.249 





0.25 


10 


10 


7 


6 


0.381 


0.014 


0.66 


0.026 





0.25 





0.25 


10 


15 


8 


7 


0.389 


0.012 


0.745 


0.06 





0.25 





0.25 


10 


20 


9 


10 


0.0] 2 


0.001102 


0.3G9 


0.017 





0.25 





0.25 


10 


25 


10 


20 


0.363 


0.019 


0.06 


0.194 





0.25 





0.25 


15 


5 


6 


5 


0.439 


0.003772 


0.728 


0.052 





0.25 





0.25 


15 


10 


7 


6 


0.723 


0.05 


0.715 


0.046 





0.25 





0.25 


15 


15 


8 


7 


0.535 


0.001225 


0.608 


0.012 





0.25 





0.25 


15 


20 


9 


10 


0.66 


0.025 


0.628 


0.016 





0.25 





0.25 


15 


25 


10 


20 


0.639 


0.019 


0.178 


0.104 





0.25 





0.25 


20 


5 


15 


4 


0.505 


0.00002566 


0.406 


0.008777 





0.25 





0.25 


20 


10 


16 


9 


0.541 


0.001673 


0.609 


0.012 





0.25 





0.25 


20 


15 


17 


14 


0.472 


0.0007671 


0.00226 


0.248 





0.25 





0.25 


25 


10 


20 


9 


0.487 


0.0001627 


0.551 


0.002556 





0.25 





0.25 


25 


15 


22 


14 


0.509 


0.00008313 


0.487 


0.0001641 





0.25 





0.25 


25 


20 


24 


19 


0.528 


0.0007997 


0.576 


0.00584 





0.25 





0.25 


50 


25 


-1 


24 


0.791 


0.085 


1.68 


1.392 





0.25 





0.25 


50 


30 


6 


29 


0.752 


0.063 


1.144 


0.415 





0.25 





0.25 


50 


40 


9 


39 


0.755 


0.065 


0.023 


0.227 





0.25 





0.25 



a = 7,P = TandR = 0.5 



Where both a and (3 are the scale parameters for X and Y respectively. R is the 
reliability function ,also. MSEl is The mean square error of Ri , MSE2 is The mean 
square error of R2 , MSE3 is The mean square error of R^ and MSE4 is The mean square 
error of i?4 . 
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